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We present a new approach to the analysis of quark-binding effects in inclusive decays of heavy 
mesons within the relativistic dispersion quark model. Various differential distributions, such as 
electron energy spectrum, q 2 - and Mx -distributions, are calculated in terms of the B meson soft wave 
function which also determines long-distance effects in exclusive transition form factors. Using the 
quark-model parameters and the B meson wave function previously determined from the description 
of the exclusive b —* u transitions within the same dispersion approach, we provide numerical results 
on various distributions in the inclusive B — > X c £i>e decays. 



I. INTRODUCTION 



' Inclusive decays provide a promising possibility to determine the CKM matrix elements describing the mixing 
of b quark, since a rigorous theoretical treatment of these decays, including nonperturbative effects, is possible. A 
consideration based on the Operator Product Expansion (OPE) and the Heavy Quark (HQ) expansion jjj allows one 
to connect the rate of the inclusive B meson decay with the rate of the b quark decay. An important consequence 
of the analysis based on the OPE is the appearance of the HQ binding effects in the integrated rates, both total 
and semileptonic (SL), of heavy meson decays only in the second order of the HQ expansion g]. These second order 
I . corrections are expressed in terms of the two hadronic parameters, Ai and A2. The latter are the mesonic matrix 
' elements of the operators of dimension 5 which appear in the OPE of the product of the two weak currents. The 
differential distributions can be calculated as expansions in inverse powers of HQ mass txiq . 

Whereas presumably providing quite reliable results for the integrated SL decay rate, the OPE method encounters 
difficulties in calculating various differential distributions. For instance, before comparing the OPE-based results for 
the differential distributions with the true distributions a proper smearing over duality interval is necessary. 

There are several reasons which yield complications in the calculation of some differential distributions, arising 
. mostly in the resonance region near zero recoil, namely: 

1. the duality-violating l/m,Q effects (i.e. the difference between the true distributions and the smeared OPE results) 
Oh, in the differential distributions near zero recoil which originate from the delay in opening different hadronic channels, 
■ as noticed in ||. Although these effects are cancelled in the integrated SL rate, they can considerably influence the 
' kinematical distributions near the zero recoil point; 

■ 

• • . 2. the convergence of the OPE series for the differential distributions persists only in the region where the quark 
produced in the SL decay is sufficiently fast. This means that the OPE cannot directly predict distributions in some 
, kinematical regions, such as: 
;_i ' 

• the photon energy spectrum dT/dE-y in the radiative B — > X s ^: the window in the photon energy between 
mq/2 and Mq/2 turns out to be completely inaccessible within the OPE formalism [H; 

• the lepton energy spectrum dT/dEi at large values of Eg in semileptonic or rare leptonic decays; 

• the lepton (^-distributions in SL B — > X c (X u ) and rare B — > X s decays at large q 2 near zero recoil; in this 
region one encounters both the quark-binding and duality-violating effects. 

Problems related to the quark-binding effects can be solved in principle by performing proper resummation of 
the nonperturbative corrections which in practice however leads to the appearance of a priori unknown distribution 
functions ||[r|. 

The inclusion of the quark-binding effects in heavy meson decays was first done in Q , where an unknown distribution 
function of a heavy quark inside the heavy meson was introduced. Evidently, this distribution function is connected 
with the wave function of the heavy meson which also determines the exclusive transition form factors. To put this 
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connection on a more solid basis, it is reasonable to consider the inclusive and exclusive processes within the same 
approach. 

We argue in this paper that the constituent quark model (QM) can be used an efficient tool for calculating differential 
distributions in inclusive decays of heavy mesons, covering also kincmatical regions where OPE cannot provide a 
rigorous treatment. Namely, the constituent quark model allows one to take into account quark-binding effects 
in inclusive heavy meson decays in terms of the meson soft wave function. The latter describes the heavy meson 
properties both in exclusive and inclusive processes and thus allows one to consider on the same ground long-distance 
effects in various kinds of hadronic processes. 

Quark-model calculations of inclusive distributions are essentially based on the evaluation of the box diagram (see 
Fig. la later on) by introducing the heavy meson wave function in one way or another. To illustrate the basic features 
of such an approach as well as its advantages and limitations it suffices to consider the case of a nonrelativistic (NR) 
potential model with scalar currents. Inclusion of relativistic effects can be then performed. 

Let us consider a weak transition induced by the scalar current J = cb, where both b and c are heavy. To make the 
nonrelativistic treatment consistent we assume that 

TUb, m c 3> 5m = m b — m c 3> A, (1) 

where A is the typical scale of the quark binding effects in the heavy meson. In the NR theory the general expression 
for the hadronic tensor 

W{q ,q) = ~Im J < B\T(J(x)J+ (0))\B > e~ lqx dx (2) 

is reduced to the form 

W(q , q) = -7m < B\G cd (M B - q° - iO, q)\B > . (3) 

7T 

Here G cd (E, q) = (H cd (q) — E)^ 1 is the full Green function corresponding to the the full Hamiltonian operator of the 
cd system with the total momentum q 

* C M=rn c + md + ^ + ^ + V cd (f). (4) 

Thus, the hadronic tensor is the average of the full cd Green function over the ground state of the full bd Hamiltonian 

H bd \B >= M B \B >= (m b + m d + e B )\B > . (5) 
In the rest frame of the B-meson one has 

k 2 k 2 

H bd = m b + m d + h h V bSi {r) = m b + m d + h bd (6) 

lm b lm d 

The following relation provides a basis for performing the OPE in the NR potential model M 



^2 



H cd --(M B - q ») = -(8m-^--q°) + h bd -e B + I - ) - + O I ^ ) (7) 



2m r 



k 2 + V\ ( 1 1 \ kq ( 7 3 <5m 



ra c m b 



where 7 ~ Aqcd and we have assumed the following expansion of the potential 

Starting with (^) one constructs an OPE series using the amplitude of the free b — > c quark transition as a zero-order 

approximation (hereafter referred to as the standard OPE). By virtue of the equations of motion, (h bd — e B ^ \B >= 

one observes the absence of 1/toq corrections to the leading order (LO) b — > c amplitude, so that the corrections 
emerge only at the I/toq order. Being completely reliable for the calculation of the integrated decay rate, this choice 
of the zero-order approximation turns out to be inconvenient however for calculating differential distributions. In 
particular, the distribution in the invariant mass of the produced hadronic system, Mx, becomes very singular and 
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is represented via S(Alx — m c ) and its derivatives, such that the 1/toq corrections are even more singular than the 
LO result. This is the price one pays for the choice of the zero-order term. 

It is clear that the free-quark transition amplitude is not the unique choice of the zero-order approximation, at least 
in quantum mechanics. For instance, another structure of the expansion can be obtained if the free cd Green function 
is used as the zero-order approximation. 

In the NR quantum mechanics the relation between the full and the free Green functions is well known and reads 

G-\E)=H-E, Go 1 (E) = H -E, G~ l (E) — Gq 1 (E) = V, (9) 

or, equivalently, 

G(E) = G (E) - G Q (E)VG(E) = G (E) - G (E)VG (E) + G Q (E)VG Q (E)VG Q (E) + ... (10) 

For the heavy quark decay, in most of the kinematical g 2 -region except for a vicinity of the zero recoil point, the 
Green function Go behaves as 1/toq, and, since the matrix elements of the operator V remain finite as toq — > oo, 
the series ([lO]) is an expansion in powers of 1/rriQ. Notice that in the NR potential model the expansion (|l0|) is fully 
equivalent to the OPE series obtained from Inserting the expansion ( jToj ) into the expression for the hadronic 
tensor W given by eq. (ft), we come to the series shown in Fig 1. 
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FIG. 1. Expansion of the hadronic tensor in the quark model: (a) - the box diagram which provides the LO free-quark 
contribution, (b,c) - diagrams contributing in subleading I/tuq orders and containing the final state interaction V = V c j. The 
short-dashed lines represent W bosons. 



The LO term is the box diagram of Fig. la with the free c and d quarks in the intermediate state. The corresponding 
analytical expression reads 

W QM = Il m < B\G° cS (M B - q°, q)\B > . (11) 

This is the quantity usually taken into account in QM calculations. It is easy to see, that the SL decay rate calculation 
based on the box diagram of Fig. 1 only, reproduces the free quark SL decay rate in the HQ limit, but should contain 
also 1 / ttiq corrections (see also the general structure of the QM results of Ref . |l(J ) . Namely, the next-to- leading 
order (NLO) term in the expansion ( |l0| ) is the diagram with a single insertion of the potential between the free c and 
d quarks (the diagram of Fig. lb). It has the order 1/mq and precisely cancels the 1/toq contribution of the QM box 
diagram, yielding in this way the absence of the 1 /toq correction in the difference between the decay rates of bound 
and free heavy quarks. 

Thus, the QM box-diagram calculation is just the first-order term in an alternative expansion of the full Green 
function: unlike the standard OPE series which starts from a single c-quark in the intermediate state, the QM starts 
from the free cd pair which is the eigenstate of the Hamiltonian 

„ >1 = mt+mj + <^: + |_. (12) 

Hereafter we refer to the expansion of the decay rate based on the expansion ( |l0| ) of the Green function as the QM 
expansion. 

Summing up, the quark model provides an alternative 1/toq expansion with the following properties: 

1. the box diagram of Fig. la provides the LO 1/toq term and reproduces the free-quark decay in the limit 
toq — ► oo. All other terms contribute only in subleading 1/toq orders; 
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2. the differential distributions in any I/tuq order are convergent for almost all allowed q 2 , except for the region 
close to zero recoil; 

3. before comparing the calculated differential distributions based on the expansion ( |l0| ) with the true distributions 
in the resonance region a proper smearing over some duality interval is required; 

4. the 1 / rriQ correction to the LO term is nonvanishing. 

Clearly, the properties 1-3 are completely equivalent to the standard OPE, while the property 4 makes the QM 
expansion much less convenient than the standard OPE, at least for the calculation of the SL decay rate. However: 

5. the expansion ( |To| ) turns out to be more suitable for the calculation of the differential distributions, e.g. for the 
calculation of dT/dMx- in this case the LO result (the box diagram of Fig. la) is well-defined in the whole 
kinematical region as well as the higher order corrections to it. Thus already the box diagram is appropriate for 
comparison with the experimental dT/dMx at all Mx apart from the resonance region. Beyond the resonance 
region no additional smearing of the calculated dT/dMx is required. 

In full QCD the situation is of course much more complicated and hadrons are coherent states of infinite number 
of quarks and gluons. Nevertheless, many applications of the constituent quark model have proved the treatment of 
mesons as bound states of two constituent quarks to provide a reasonable description of their properties. From this 
viewpoint the arguments given above remain valid. Namely, the box diagram represents the main contribution to 
the hadronic tensor which reproduces the free-quark decay in the infinite quark mass limit. However hadronic tensor 
calculated from the box diagram contains the 1/tuq term compared with the free-quark decay tensor. This linear 
1 /m,Q term is known to be cancelled by the 1 / wiq order contributions of higher order diagram and to be absent in 
the full expression. In practice, however the 1/toq term of the box diagram is not so dangerous: namely, the hadronic 
tensor calculated from the box diagram, as well as all corrections given by the other graphs, are regular in the whole 
kinematical region. Thus the box diagram should provide a reasonable description already appropriate for comparison 
with experiment. Moreover, the box-diagram result can be further improved by effectively taking into account the 
higher order term which kills the 1/mg correction contained in the box diagram. We follow this strategy in our 
analysis and perform a relativistic treatment of quark-binding effects within a constituent quark picture. 

Our consideration of the quark binding effects in inclusive SL decays is based on the relativistic dispersion formu- 
lation of the quark model previously developed for the description of meson transition form factors |l l| . Within this 
approach, the inclusive decay rates as well as the exclusive hadron transition form factors are given by double spectral 
representations in terms of the soft meson wave functions. The double spectral densities of these spectral represen- 
tations are obtained from the corresponding Feynman graphs. The subtraction terms in spectral representations for 
exclusive transition form factors are fixed by requiring the structure of the HQ expansion in the QM to match the 
structure of the HQ expansion in QCD. In this paper we proceed along the same lines in inclusive processes. 

Our main results are the following: 

• we construct the double spectral representation of the hadronic tensor within the constituent quark model 
starting with q 2 < 0. The hadronic tensor is represented in terms of the soft wave function of the B meson and 
the double spectral density of the box diagram. The hadronic tensor at q 2 > is obtained by the analytical 
continuation. Then the I/tuq expansion of the spectral representation of the decay rate is performed and the 
LO term is shown to reproduce the free quark decay rate. The subtraction is defined in such a way that the 
1/roq correction to the SL decay rate is absent. This corresponds to effectively taking into account other terms 
beyond the box-diagram approximation which contribute in subleading 1/mg orders. Moreover an approximate 
account of the 1/rrig effects of the whole series within the box-diagram expression is possible. This is done by 
introducing a phcnomcnological cut in the double spectral representation of the box-diagram which affects only 
the differential distributions at large q 2 . This cut brings the size of the 1/toq corrections in the T(B — > XJLvi) 
in full agreement with the OPE result and keeps the LO and 1/rriQ correction unchanged. The cut yields 
differential distributions which are finite also in the endpoint g 2 -region where the HQ expansion series is not 
properly convergent; 

• various differential distributions are calculated in terms of the £?-meson soft wave function. These distributions 
are regular in the whole kinematically accessible region and, apart from the resonance region (where the exact 
distributions are dominated by single resonances and a proper smearing over the duality intervals is necessary), 
can be directly compared with the observable values. The main effect of quark-binding upon these distributions is 
determined unambiguously through the soft wave function, while the \jm 2 Q corrections depend on the particular 
details of an approximate account of the higher-order terms in the series (|l0|). However, in practice these details 
are not essential due to the following two reasons: first they are numerically small, and second, the size of the 



4 



1/rriQ corrections in the integrated SL rate is close to the OPE result. So we expect that the size of the 1/rng 
corrections is reasonably reproduced also in other quantities; 

• we perform numerical estimates of various differential distributions in inclusive decays in terms of the B-meson 
wave function known from the description of the exclusive processes [ ^2| , p^[ . 

The paper is organized as follows: in the next section we present necessary formulas for the free-quark decay and 
also the OPE prediction for the total integrated rate up to 1 /rrig corrections. In Section 3 we construct the dispersion 
representation for the box diagram at q 2 < and discuss its analytical continuation to q 2 > 0. Section 4 performs the 
1/mq expansion of the hadronic tensor in the quark model, and Section 5 presents numerical results for differential 
distributions. Finally, a brief summary is given in the Conclusion. 



II. FREE QUARK DECAY AND OPE 

Let the effective Hamiltonian governing the quark transition Q — ^ Q' with the emission of the particle <fi have the 
following structure 

H eff {x) = Q'{x)YQ{x)(j){x) (13) 

where T denotes a relevant combination of the Dirac matrices. Following notations of ref JlTf for the exclusive form 
factors, we denote the parent heavy quark Q also as Qi and the daughter quark Q' as Q\. 

A tree-level rate of the free-quark decay initiated by this effective Hamiltonian averaged over the polarizations of 
the initial quark Q 2 and summed over polarizations of the final quark Q± has the form 

dTo (2*)« f ]n? dk0 m _ _ q2)s{k2 _ ki _ ^ 



dq 2 2rri2 

(2ir) 4 nX^imlmlq 2 ) 



?o| 2 , (14) 



2m 2 (27r) 6 2m\ 
where q 2 is the mass squared of the particle <f> and 

|T | 2 = \ ^( fc 2)f ( Wl + ki)£u a {k 2 ) = Jsp ((m 2 + fe 2 )f (m x + fc x )f ) . (15) 



Hereafter we use the notation X(x, y, z) = {x + y — z) 2 — Axy. 

These formulas can be readily applied to the particular cases of radiative and SL decays. In the latter case one 
needs to multiply dT^/dq 2 by the the leptonic tensor L(q 2 ) to obtain the full differential distribution dTQ L /dq 2 = 
dTg/dq 2 ■ L(q 2 ). Hereafter the inclusion of the leptonic tensor in the definition of dTo/dq 2 is understood and we drop 
the superscript SL. 

In the case of the SL b — > clvi transition the free-quark differential decay rate reads explicitly as (cf., e.g., |§J|) 

^0 G^IVbcl 2 1 , 1/2/ 2 2 2\ril 2 2 2\ 

~dq* = 96tt 3 rr? ^ m ^ m c^ ) c ( m b, m c:<l ). (!6) 
where Gf is the universal Fermi constant, Vb c is the CKM matrix element, and 

C(ml m 2 c , q 2 ) = (m 2 - m 2 c f + q 2 (m 2 + m 2 c ) ~ 2q 4 . (17) 
The integrated SL rate is then given by 

r(m b -m c f ^ G 2 ]Vbc] 2 
T0 = J dq d^ = ^^' mb ' l0ir) (18) 

with Io(r) = 1 — 8r + 87- 3 — r 4 — 12r 2 ln(r), r = (mc/mi,) 2 . 

The OPE predicts that in the decay of the heavy meson, the l/m^ corrections to the free quark rate ( |l8| ) are absent, 
and up to 1/m 2 terms the integrated SL rate is given by (cf., e.g., 0,0]) 



To 



Ai+3A 2 A 2 (1-r) 



1+ 2mg 6 m 2 Jo(r) 



(19) 
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Here Ai and A2 are the hadronic matrix elements of the operators of dimension 5 appearing in the OPE of the product 
of the two weak currents. The value A2 = 0.12 GeV 2 is well known from the B — B* mass splitting, whereas the 
knowledge of Ai is loose and present estimates range from —0.6 GeV 2 to 0. Note that in the NR quark potential 
model one has Ai = — < k 2 >, where k is the relative momentum of the constituent Qq pair (cf, e.g., Typically, 
the NR quark model estimates of Ai range from —0.6 to —0.4 GeV 2 . 



III. INCLUSIVE MESON DECAY IN THE QUARK MODEL 



We now proceed to the calculation of the inclusive rate for the decay of a pseudoscalar meson with mass Mi 
containing a heavy quark Q 2 , which we will refer to as Pq 2 , induced by the quark transition dl3). 

We start with the box diagram of Fig. la. Our notations shown in Fig || follow those of rcf where transition 
form factors within the similar dispersion approach have been considered. 
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FIG. 2. Momentum notations in the box diagram. 



dT(P Q2 -> X Ql <t>) _ (2tt) 4 f dk 1( lk 3 dk^ 



The decay rate corresponding to this diagram can be written in the form 

, M/( ; ( , V)I) ^T\ 2 5(k 2 -m 2 )S(k 2 - m 2 )5(k 2 -q 2 )S(p 1 -k 2 -k 3 -k 4> ), (20) 

where T is the amplitude of the meson decay Pq 2 — ► Qiq 3 4>. 

In the dispersion approach the decay rate can be written in the form of the following spectral representation: 

dT(P Q2 -» XqJ) (2^) 4 f d Sl ds 2 ir\ 1 / 2 (si,S2, q 2 ) ~ AI 2 , _ 
dq~ 2 ( ai -M?)2 2^ A (si,S2,q) (21) 

where the spectral density A(si, s 2 , q 2 ) is connected with the double discontinuity of the box diagram shown in Fig. 
1. 

The amplitude A of Fig. la depends on six Lorentz scalar variables p 2 ,p' 2 ,p 2 , q 2 , q' 2 , n 2 , where the momenta satisfy 
the following relations p\ — p[ = q — q' = k, and p\ — q = p\ — q' = p 2 - The amplitude can be written as double 
spectral representation in p\ and p' 2 

A(p\,p' 2 ,pla 2 ,q' 2 ,K 2 ) = J s * S J s ^ 1 p? A(s 1 , s ' 11 plq 2 ,q' 2 , K 2 ). (22) 

In this equation A is the full spectral density of the amplitude and includes properly defined subtraction terms. The 
spectral density A is calculated from the double discontinuity Ad, 

Ad = 1 disc^discy ,4(s;l, s'i,p|, q 2 ,q' 2 ,n 2 ), (23) 

through a relevant subtraction procedure. 
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A. The spacelike region 



The double discontinuity Ajj of the box diagram of Fig. la can be calculated at q 2 < and q 12 < by placing all 
intermediate quarks on their mass shells but keeping the initial and final mesons virtual and having the masses squared 
si and s[, respectively. To obtain the corresponding expression at positive values of q 2 = q' 2 which is necessary e.g. 
for SL decays, we shall perform the analytical continuation in q 2 . At q 2 , q' 2 < the procedure explained in detail in 
fj"lf yields the expression 

Ad(si, s'i, s 2 , q 2 , q' 2 , n 2 ) = (n x q f dk\dk2dk' 2 dkz5{k 2 — m\)5(k 2 — m 2 )S(k 2 2 — m 2 )S(k 2 — to 2 ) 



(2tt) 9 

x5(pi - k 2 - k 3 )S(p[ - k! 2 - k 3 )S(p2 - ki - k 3 ) 



x(-l)Sp (i 75 G{ai) (m 3 - k)h 5 G[s[) (m 2 + fc 2 )f (mi + fc x )f (m 2 + k 2 )) , (24) 



where the momenta satisfy the following relations p[ = p± + k, q' = q + K, p~2 = Pi — q, V\ = s ii Pi — s 'ii Q 2 = 9 2 i 
q' 2 = q' 2 , p 2 = S2 ■ Notice that the quark structure of the (virtual) pseudoscalar meson transition into two real quarks 
is described by the vertex G(si)Q2(£;2)*75<Z3(&3)<5(pi — fea — k 3 ), where p 2 = si, k 2 = m 2 ,, fc| = m§. 

The spectral density A{s\, S2, q 2 ) in (|2^) is connected with the forward double spectral density of the box diagram 
A(si,s' 1 ,S2,q 2 ,q' 2 ,K 2 ) as follows 

lim A(si, s' 1; 82, q 2 , q 2 , k 2 ) — S(si - s[)A(si, s 2 , q 2 ). (25) 

K 2 — *0 

The form of the subtraction procedure in the spectral representation (^) cannot be determined within the dispersion 
approach and should be fixed from some other arguments. We can determine the subtraction term by requiring the 
absence of the l/mq corrections in the ratio of the bound and free quark decay rates. As we have discussed, this 
corresponds to taking into account the 1 /tuq terms of other diagrams which are known to cancel the 1 /mq term in 
the bix diagram. Our subtraction prescription explicitly reads 

A(s 1: S2,q 2 ) = -±A D (8 1 ,8 2i q 2 ) (26) 



where Ap is connected with the double discontinuity of the forward amplitude 

lim A D (s 1 ,s' l ,S2 1 q 2 ,q 2 ,K 2 ) = 5(s\ - s' l )A D {si, s 2l q 2 )- (27) 

Explicit calculations give for Ad the following expression 

^( Sl , S2 , g 2 ) = 2Al/2 ^ i!g2) (-S P ). (28) 

Notice that the argument of the 9 function in eq (p8| ) is just the same as for the spectral density of the triangle graph 
describing the form factor in the dispersion approach which can be read off from [pi. Since all quarks are on their 
mass shell, the trace can be rewritten in the following form 

-Sp = 2|T | 2 ( Sl -(m 2 -m 3 ) 2 ), (29) 

where |Tb| 2 is just the square of the free-quark amplitude of eq jl5|). Finally, isolating the free-quark decay amplitude 
and solving the kinematical function we come to the following dispersion representation for the differential inclusive 
decay rate at q 2 < 



rfr ts ( 2\ f j g2 ( s i) si-(m 2 -TO 3 ) 2 m 2 



i+(si 
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-~K {q ) / d Sl - —2-2 —J 7= u727 2 2-^T / ds ? X S 1' S 2,9 , (30) 

dq 2 J {si-M(Y 8tt j si ^/si A V 2 [mi, mi, q 2 ) J s -{s u q 2 ) 

where 

K( 2, 1 dT Q _ G 2 F \V 2l \ 2 C{mlm 2 ,q 2 ) 

° Kq } ~ W 2 [mlmlq 2 ) dq 2 96tt 3 m\ ' K ' 
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The limits s 2 are obtained by setting 77 = ±1 in the equation 



sa(si, q 2 ) = (mi + m 3 ) 2 + — (si - (m 2 + 7713) 2 ) + — (ui - l)(si -m\- ml) 

m 2 ' m 2 

+ r/ — -A 1 / 2 (si, m|, m 2 )^Ju 2 — 1, 
m 2 

where the quark recoil w is defined as follows 

q 2 = (m 2 — mi) 2 — 2?tt,ito 2 (cl; — 1). 



Tfl\ 



(32) 



(33) 



Note that in eq (30) the free-quark differential rate dT^/dq 2 factorizes out, so that the differential rate for a bound 
quark is a product of the free-quark differential rate and a bound state factor, as already noted in |LQ |. Hereafter we 
use the notation ip(s) = G(s)/(s — M 2 ). The normalization condition of the soft wave function ip(s) obtained from 
the elastic vector form factor of the heavy meson at q 2 — reads pT| 



dsi<p (si) Isi - (m 2 - m 3 ) — -= = 1. 

L J 87T Z Sl 



(34) 



It is convenient to rearrange eq (30) by isolating under the integral the structure similar to the structure of the 
normalization condition (p4[) 



= K Q (q 2 ) J ds 1 ip 2 (s 1 ) [si - (m 2 - m 3 ) 2 ] 



AV 2 ( ai ,ml,m§; 
8tt 2 si 



p{si,q 2 ), 



p{si,q 2 



m 2 



si A 1 / 2 (si, m|, m|)A 1 / 2 (m2, m 2 , g 2 ) 



ds 2 X^ 2 (s u s 2 ,q 2 ). 



(35) 



As we shall see later, p(si,q 2 ) ~ A 1 / 2 (m 2 , mf , g 2 ) in the HQ limit, so that thanks to the normalization ([mJ) of the 
soft wave function one gets dT/dq 2 — > dTo/dq 2 as mQ — > 00. 



B. The timelike region and the anomalous contribution 

To obtain the spectral representation at q 2 > we perform the analytical continuation in q 2 . This procedure is 
done along the same lines as in the case of the transition form factor which has been discussed in detail in . As a 
result of this procedure in addition to the normal part which is just the expression (35) taken at q 2 > 0, an anomalous 
part emerges due to the non-Landau type singularities of the Feynman graph. Thus, in the region q 2 < (m 2 — mi) 2 
the representation for p(si 1 q 2 ) takes the form 



m 2 



iTA 1 / 2 ( Sl 



,m 2 ,m3 



)A!/ 2 ( r 



f 2 . - 



lq 2 ) 



s+(si,9 ) 



s 2 ( fi ii9 2 ) 



ds 2 \ 1 / 2 (s 1 , S2 ,q 2 ) + 26(q 2 )9(s 1 > fl ?) 



s 2 (ai,q ) 



ds 2 X 1/2 (si,s 2l q 2 ) 



s 2 ( S 1.9 2 ) 



(36) 



where 



S 2 (si,q 2 ) = {Vs 



^f^?) = 



q 2 ) 



q + m 2 - m\ 
2 V? 



\ 



q 2 + m\ — m\ 



+ m\ — mk. 



(37) 



The g 2 -behavior of the anomalous term is determined by the lower limit of the s\ integration, s\(q 2 ). Namely, its 
contribution to the SL rate reads 



I (£Y^ anorn 

r du> 



A 3 



mij^/uj — 1 



R° 



(38) 



where 
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R anom (uj) = J dkH 2 B {k 2 ). (39) 

^i-(w-i) 

A 2 v ; 

Here \k\ = A 1/,2 (si, m 2 , m 2 )/2y/si is the relative momentum of the quarks inside the B-meson, and the wave function 
\& s(fc 2 ) can be expressed through ip(s) jll). In terms of ^#(fc 2 ) the normalization condition ( |34| ) takes the form 
/ ^ 2 B (k 2 )dk 2 — 1 such that R anom (l) = 1. Since the soft wave function is steeply falling beyond the confinement 
region where k 2 < A 2 , the anomalous contribution becomes inessential already at ui — 1 ~ A/ mi,. Only in the 
endpoint region lo — 1 < the anomalous contribution to the differential distribution becomes strong and diverges 

m Q 

ilKe r — ^ ^/ztt • 

The contribution of the anomalous term to the integrated SL rate is Y anom /Y ~ A 2 /toq and comes from the 
endpoint region, whereas the rest of the phase space < q 2 < A/tuq provides only the relative A 3 /toq anomalous 
contribution to the SL decay rate. 

A 2 

Therefore, the anomalous contribution is negligible at all to except for the endpoint region oj — 1 < -V, which is m 

fact a very narrow region near zero recoil. As we have discussed, the HQ expansion for the differential distributions 
is anyway ill-defined in this kinematical region. Contributions of the same order of magnitude come also from other 
terms in the expansion ([lO]), and keeping this anomalous contribution is beyond the accuracy of our considerations. 
Thus we shall systematically omit the anomalous contribution in numerical calculations. 



IV. THE HEAVY QUARK EXPANSION OF THE INCLUSIVE DECAY RATE IN THE QUARK MODEL 

In this section we perform the HQ expansion of the meson inclusive decay rate. We show that: 

a. in the LO the heavy meson inclusive decay rate is equal to the free quark decay rate; 

b. our subtraction prescription (|2^) leads to the differential distribution dY/dq 2 given by eq ([35]) which satisfies the 
relation (dT(B — > X c £D£)/dq 2 )/(dY(b — ► c£Di)/dq 2 ) = 1 + 0(l/m|>) in most of the q 2 region except for a close 
vicinity of zero recoil point. This property guarantees the absence of the 1/toq in the ratio of the integrated 
rates, i.e. T(B -> X c lui)/Y{b -> clvi) = 1 + 0(l/m 2 Q ); 

c. the size of the I/tuq corrections can be tuned such that they become numerically close to the OPE prediction. 
This is done by introducing the cut in the spectral representation of the decay rate of the B meson. This 
cut affects only the differential distribution dT(B — > X c lvi)j ' dq 2 at large q 2 near zero recoil, i.e. in the region 
w < 1 + 0(l/m Q ). 

The most important feature of the whole approach is that already the zero order expression provides a realistic Mj- 
distribution. Modifications b) and c) while affecting the total rate and the ^-distributions at large q 2 , only moderately 
affect the Mx-distribution, so that the latter is mostly determined only by the soft B-meson wave function. 



A. Soft wave function and normalization condition 



First, we need to specify the properties of the soft meson wave function. A basic property of the soft wave function 
<p(s, mQ, mq, A) is its strong peaking in terms of the relative quark momentum in the region of the order of A ~ Aqqd- 

For elaborating the 1 / toq expansion, it is convenient to formulate such peaking in terms of the variable z such that 
s = (mq + ni3 + z) 2 (hereafter we denote the mass of the light spectator quark as m^). In the heavy meson, the 
variable z is related to the relative quark momentum as follows 

k 2 = z{z + 2m 3 )+0(l/m Q ). (40) 

Hence, a localization of the soft wave function in terms of z means that the wave function is nonzero as z < A. In 
the heavy meson case we imply that mQ 3> m$ ~ z ~ A. 

The normalization condition ( |34] ) is a consequence of the vector current conservation in the full theory and it 
provides an (infinite) chain of relations at different l/m<g orders. Namely, expanding the soft wave function in I/iuq 
as follows 
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(p(s,m Q ,m 3 ,A) 



fmQ 



-_(f) (z,m 3 ,A) 



l + £^Xi(z,m 3 , A) +0(l/m%) 



we come to the normalization condition in the form 



dz(f> 2 l (z)^(z + 2m 3 ) 3/2 



m 3 m 3 
-Xl{z) - ~ + ■ 



2toq' 



2mc 



1. 



(41) 



(42) 



This exact relation is equivalent to an infinite chain of equations in different 1/mg orders. Lowest order relations 
take the form 



/ dz 4>l{z)^{z + 2m 3 f /2 = 1, 

/ dz <pl(z)y/z(z + 2m 3 ) 3/2 xi(z) = 1, etc. 

In particular, the Isgur-Wise function is given by the following expression through 4>q 

i 

,/ x f I — t r f drj . ( 2m 3 + z\ + Z2 

£{u)= / dzi(j)o(zi)y/zi(zi + 2m 3 ) / — 0o(Z2j I m 3 -\ — — 



(43) 
(44) 



(45) 



The expression for Z2 through z\ and r\ (— 1 < r\ < 1) is obtained by expanding (32) in 1/toq 

zi = z\ + (zi + m 3 )(ui - 1) + t]^/zi(zi + 2m 3 )y / ui 2 - 1 + 0(1/to q ), 
and for the calculation of the IW function only the LO part of this relation should be used. 



(46) 



B. The HQ expansion 

First let us consider the HQ expansion of the unsubtracted quantity pp. The normal part of pd(si, q 2 ) reads 

- 1 *?u. 



2 m 2 f drj 1/2 2 



(47) 



This representation is a convenient starting point for performing the HQ expansion. Notice that although the inte- 
gration in rj can be easily performed, it is more convenient to work out the HQ expansion before the integration. 

Assuming that m 2 is large and that the meson wave function is localized in the region Z\ ~ A we obtain the 
following expression for A(si, S2, q 2 ) valid at all q 2 



\l 2n 4 (1 , z + m 3 
\(s 1 ,s 2 ,q ) — ► m 2 H 

V m 2 

z(z + 2m 3 ) 



2 r 



A(l, «f , f 2 ) + %(z)V^ + 2m 3 )(l + <f - p 2 )A 1 / 2 (l, g 2 , f 2 ) 
m 2 



+ 



(1 + q 2 - f 2 ) 2 + ^z(z + 2m 3 )A(l, q 2 , f 2 ) 



(48) 



where g 2 = q 2 /m 2 , f — mx/mz and x( z ) = 1 — ( z + Tn 3 )/2m2- In the limit m 2 — > 00 we can expand the A(si, S2, g 2 ) 
in powers of l/m^. Notice however that an actual expansion parameter is not l/m-2 but rather 



y/z(z + 2m 3 ) 
m2A 1 / 2 (l, q 2 , f 2 ) 



(49) 



and the averaging over the B meson state implies yz(z + 2m 3 ) ~ A. Hence the region where the expansion is fastly 

converging is m^A 1 ' 2 ^, g 2 , f 2 ) ^> A. This relation can be written as \k\\ = A 1 / 2 (m 2 ,,rn 2 ,q 2 )/2m2 A, which means 
that in the rest frame of the b quark the daughter quark has a 3-momentum much bigger than A. 
The final expression reads 



p D ( Sl ,q 2 ) - AV3( m ?,mI,<7 2 )^ (1 + 

Mi V m 2 



z(z + 2m 3 ) 



2ml 



8q 2 



3A(l,g 2 ,f 2 



(50) 
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The 1/toq term in the ratio of the bound to free quark distributions is generated by the (rri2 + z + m$)/Mi term in 
PD- 
As we have discussed this linear 1 /mo term contained in the box diagram cancels against the 1 /tuq terms coming 
from other terms in the expansion (0). Thus the main contribution of these other terms in the series (ficj ) can be 
taken into account by performing the subtraction in the spectral representation of the box diagram which kills the 
1/niQ term as follows: 

p(si,q 2 ) = ^pn(si,q 2 ). (51) 



In the HQ limit z(z + 2ms) = k 2 , so after performing the subtraction we come to the following relation 



, _ dT(B -> X c iv t )/dq 2 



R( - q ) - dT (b -» cfi/^/d^ 1 + 2m| V + 3A(1, q 2 , f 2 ) ) ' (52) 

This expansion is valid in the region of q 2 such that A(l, q 2 ,f 2 ) — 0(1), i.e. in most of the q 2 phase space except for 
the region near zero recoil where A(l, q 2 , r 2 ) ~ 0. 
The expression (pi) has the following features: 



1. in the LO the ratio R{q 2 ) is equal to one and thus the decay rate of the free and the bound quark coincide in the 
HQ limit at all q 2 . Moreover, the differential distribution also coincide within the 1/mQ accuracy in most of the 
q 2 phase space, except for the region near zero recoil. This guarantees the absence of the 1/mQ corrections in 
the ratio of the integrated rates as well. Thus, our description is in full agreement with the OPE results within 
the 1/mQ order; 

2. since the box diagram represents only a part of the 1/toq corrections, we cannot expect the box diagram alone 
to reproduced correctly the l/rUg term in the ratio of the integrated rates T/T - In fact, the sign of the 1/toq 
correction in eq. (|52| ) turns out to be opposite to the OPE result (cf., e.g., with the results of refs. [|||] at 
q 2 = 0). Moreover, the whole l/m 2 Q effect in the box diagram is expressed only in terms of < k 2 >, whereas 
the 1 /ttVq corrections of the whole series contains also <V%> Q] , where V± is the 1 /tuq term appearing in the 
expansion (@) of the effective potential (e.g., the chromomagnetic operator in QCD). 

We argue however that it is possible to further modify the spectral representation of the box diagram to bring the 
size of the 1/mg term developed by this modified representation in agreement with the OPE result. This procedure 
corresponds to phenomenologically taking into account the contribution of other 1 / nig terms of the expansion (|T 

Omitting the anomalous contribution as discussed previously, the differential decay rate reads 

— =K (q') da-np (si) =-3 A {si,m 2 ,m z )—— / —A (si, s 2l q )■ (53) 



where S2 depends on r\ through eq (p2|). 

Our goal is to modify the expression (53) in such a way that the LO result and the 1/ttlq correction in the 
integrated rate remain intact whereas the l/rnq term numerically reproduces the OPE estimate. We can allow a 
strong deformation of the differential (^-distribution at large q 2 near zero recoil, where the HQ expansion is anyway 
ill-defined. We can also require the 1 /tUq correction in the differential decay rate at q 2 = to exactly reproduce the 
OPE result. 

Most easily this program may be implemented through the following two steps: first, by introducing the factor 
F(s\) = 1/(1 + fc 2 /mg) which sets the 1/toq term in the differential rate at q 2 — and, second, by changing the 
upper limit in the s 2 integration in ( |53| ) to some s^^iq 2 ) which tunes the size of the 1/toq effects in the integrated 
rate: 

^=K Q (q 2 ) I d Sl <p*( Sl ) Sl ~ ~ w 3) a xV\s x ,mlrnl)^nsi) f f A 1 ^, « 3> f)9 (s 2 < s^tf)) ■ (54) 
dq z J %it z si ^/sT J -1 2 

In order not to affect the integrated rate in the LO and the 1/mQ order, s^^iq 2 ) should satisfy certain properties. 

Assume that the soft wave function 0(si) is localized in the region si < s l lax ~ [mi + ma + j) 2 where 7 is a 
constant of order A which does not scale with mQ. Let us determine q 2 , through the equation 
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4(s? a *,<l 2 )=s? ax (<l 2 ) (55) 



where sj" is the maximal value of S2 corresponding to 77 = 1 in (p2|). Furthermore, assume that s™ ax {q 2 ) decreases 
with q 2 , and take into account that s^~(si, q 2 ) is a monotonous rising function of both s% and q 2 . Then, at q 2 < q 2 for 
all si < s™ aa: one finds the relation s^(si,g 2 ) < s™ aa; (q 2 ), and thus the (^-distribution does not feel the presence of 
the cut at all. For q 2 > q 2 the cut becomes really effective and strongly influences the ^-distribution. In order these 
changes in the cut (^-differential distribution not to change the integrated rate in the LO and 1/mg order, we need 
the q 2 to be not far from zero recoil such that the corresponding ujq = 1 + 0(l/mg). Choosing the cut in the form 

\J s™ ax (q 2 ) — mi + 777,3 + a (jn2 — TUi — e — y/cf^J , (56) 

where e ~ A and a is a rising function of 7tiq, satisfies these requirements^]. 

The parameter e accounts for a mismatch between the quark and the hadron threshold, and the form of a(mo) can 
be found from fitting the size of the 1 / m,Q corrections in the integrated rate to the OPE prediction. Notice also that 
the q 2 distributions obtained through the cut expression are even more realistic than those obtained from the uncut 
spectral representations. Numerical values used for the description of the distributions are given in the following 
section. 

Mostly important for us however is that these improvements on the g 2 -differential distributions by approximate 
account of higher order graphs affect only moderately the Mx-distribution in B — » X c lvg (as well as the photon 
lineshape in the rare B — > X s ^f decay), which are thus mostly determined by the B-meson wave function. The latter 
controlls long-distance effects also in exclusive transitions. 

This property allows us to obtain a realistic energy distribution and other observables through the soft wave function 
of the heavy meson. Thus we do not need to introduce any unknown 'smearing function' describing the motion of the 
b quark inside the B meson, but rather directly calculate the effects of the b quark motion with the soft meson wave 
function. 



V. DIFFERENTIAL DISTRIBUTIONS 



In previous sections we have considered in detail the integrated rate and the differential distribution dT/dq 2 . 
We expect the obtained spectral representations of these quantities to describe some essential features of the exact 
quantities. The \jm,Q expension of these rates reproduce the known OPE results allowing to express nonpcrturbative 
parameters through the B-meson soft wave function. Interesting results can be obtained also for other differential 
distributions, such as the Mx-distribution and the lepton energy spectrum in SL decays. 

For instance, neglecting the radiative corrections in the free-quark decay, which is the LO process within the OPE 
framework, one finds the Mx distribution in the form 5(Mx — m c ). Inclusion of the \ jm,Q corrections yields a singular 
series containing derivatives of the S function. For the interpretation of these results one needs an introduction of a 
smearing function. In the quark model the Mx-spectrum is smeared because of the Fermi-motion of the b quark in 
the B meson and the shape of the Mx-spectrum is calculable through the B meson wave function. 

The meson wave function ip(s) can be written as follows JllJ 



(m? — 



V2 y/s - {m Q - m q ) 2 s-V 4 



1 One could choose a more sophisticated parameterization of S2 lax (q 2 ) to reproduce a correct q 2 -behaviour of dT/dq 2 near 
zero recoil point. For instance, taking into account that the lightest final meson is pseudoscalar, we can write y s™ ax (q 2 ) = 

mi + rns + ap(rriQ)(Mp Q — Mp , — y/q 2 ) 3 ^ 2 yielding the correct behavior near y/tf* = Mp Q — Mp q i where only one P— wave 
decay channel Pq — > Pq/Ivi is opened. In addition, in the heavy quark limit the S 1 — wave transition Pq — > VqiIvi requiring 
another functional dependence yj s™ ax (q 2 ) — mi + + av{mQ){Mp Q — Mp Ql — y/q 2 ) 1 ^ 2 is opened at yjlf- = Mp Q — Mv Q , 
with only small delay in q 2 of order A 2 . So the effects of opening this channel are even more important and should be also taken 
into account. However in the region of large q 2 with few opened channels the inclusive consideration is anyway not working 
properly, and taking into account such subtle effects is beyond the accuracy of the method. So in numerical calculations we 
proceed with the phenomenological cut provided by eq (p6|). 
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where k = \k\ = \ 1 / 2 (s,mQ,m 2 ^)/2y r s and w(k 2 ) is the ground-state S-wave radial wave function, normalized as 

As shown by the analysis of the exclusive form factors within the dispersion approach |13j , for obtaining a good 
description of the lattice results on the B — > 7T, p form factors at large q 2 it is sufficient to take into account properly 
only the confinement scale effects, i.e. to assume a simple exponential parameterization of the radial wave function in 
the form WB(k 2 ) ~ exp(— k 2 /2(3 2 3 ). Then, the quark masses m u — 0.23 GeV and m& = 4.85 GeV as well as the value 
Pb = 0.54 GeV have been determined from fitting the lattice data on the form factors. This value of (3b corresponds 
to Ai = — < k 2 >~ —0.44 GeV 2 . Finally, the mass of the c-quark is taken equal to m c = 1.35 GeV in agreement 
with current estimates of the difference m& — m c (~ 3.5 GeV). 

The parameters of the cut have been chosen according to the criteria of the previous section, obtaining e = 
(niQ — rriQ') — (Mq — Mq>) and a = 1.82 + 0.029toq, where m,Q> is the mass of the parent heavy quark and Mq* is 
the final meson lowest mass. Effectively, this means that 5^, oa . = (Mq — Mq>) 2 . The corresponding integrated decay 
rate T is plotted in Fig. |^ as a function of 1/ itiq and compared with the OPE predictions (|19|) ^| One can clearly see 
that the calculations based on both the unsubtracted pr> ( |47| ) and subtracted p ( ]5l| ) densities predict a larger rate 
for a bound heavy quark that contradicts OPE, whereas the introduction of the phenomenological cut s™ ax (q 2 ) ( |5^ ) 
brings our QM predictions in perfect agreement with the standard OPE in the whole range of considered values of 



1.05 

o 

U 

I- 1.00 
0.95 



0.00 0.05 0.10 0.15 0.20 0.25 

1 / m Q (GeV 1 ) 

FIG. 3. The ratio of the integrated rates of the bound-to- free quark SL decay Q — * Q'£i>e vs the inverse HQ mass 1/ttiq at 
fixed value of vtlqi /mq — m /mb = 0.28. Dotted line - the rate calculated from the initial spectral representation of the box 
diagram which contains I/thq correction, dashed - with the proper subtraction killing the 1/rriQ term but without tuning the 
size of the 1/rriQ effects. Solid - final result which also includes the cut bringing the size of the 1/rriQ effects in agreement with 
the OPE prediction. Upper and lower dot-dashed lines are the OPE results (hpj) corresponding to Ai = and Ai = —0.6 GeV 2 
(with A 2 = 0.12 GeV 2 ), respectively. 

Figure || shows the influence of the cut upon the q 2 -distribution for the B — * X c lT>n decay. As already discussed, the 
introduction of the (j 2 -dependent cut in the spectral representation does not change the differential q 2 distributions 
at small q 2 but it strongly affects the region of large q 2 . In particular, the cut provides the vanishing of dT/dq 2 at 
the physical threshold q 2 = (Mb — Md) 2 — 11.6 GeV 2 . Let us point out again that this cutting procedure does not 
affect the integrated rate at the leading and subleading I/ttiq orders. Note also that the differential distributions 
dT/dq 2 given by the dashed and solid lines in Fig |](a), are equal to each other at q 2 = and match the OPE result 
for dT/dq 2 (q 2 =0). 



2 To compare the calcuated integrated rate 1 as a function of tuq with the OPE result within the 1/mg order, it is reasonable 
to adopt the expansions of the hadron masses Mq and Mqi also up to the second order in 1/mg. Then one finds e = — |(Ai + 
3A 2 ) • (1/mg/ - 1/mg). Setting Ai = -0.44 GeV 2 and A 2 = 0.12 GeV 2 , yields e = 0.10/mg at m Q ,/m Q = m c /m b ~ 0.28. 
We point out that in calculations for the real B decays we use experimental values of hadronic masses (involving all orders in 
l/m Q .) 
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FIG. 4. Distribution dT /dq 2 in B —* X c £ue decays vs the squared four-momentum transfer q 2 : (a) ratio of the bound-to-free 
quark decay R{q 2 ) = (dT / dq 2 ) / (dTo / dq 2 ) , notation of lines same as in Fig. ^[ (b) Differential distribution in the bound (solid) 
and free (dot-dashed) SL quark decay. Parameters of the cut (J5(i|) are a = 1.96 and e = 0.091 GeV. 



Mostly interesting seems to be the calculated M^-distribution, reported in Fig|| Our result should be compared 
with the LO OPE result S(mx — Trie). One can see that already the box diagram of the quark model provides a smooth 
and reasonable (beyond the resonance region) distribution, which is only moderately affected by a proper account of 
the subleading 1/toq effects. At large Mx the calculated distribution does not require any additional smearing and 
can be directly compared with the experimental results. 
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FIG. 5. Distribution dT /dMx in B — > X c lT>£ decays vs the invariant mass of the produced hadronic system Mx- Lines same 
as in Fig. ^. 



The double differential distribution d 2 Y / dEdq 2 is given by the following expression 



d 2 T 



\V2i\ 2 f, 2( ^i ~ (m 2 -m 3 ) 2 u/2 , 2 2 -,H s i) f 1 d V a ( 
^T- J ds W ( Sl ) g-g- A ( Sl ,m 2 ,m 3 )-^ J [go 



> E 



dEdq 2 1287r a J "~ ir 8tt 2 si " ^^-"■■'*> m ^ J_ % 2 ~ \™ ' ' 4E 

x (2q 2 w 1 (s 1 , 8 2 , q 2 ) + [4E(q° - E) - q 2 }w 2 { Sl , s 2 , q 2 ) + 2q 2 (2E - q°)w 3 ( Sl , s 2 ,q 2 )) 



{s 2 < sT x {q 2 )) 



(58) 



In this formula q° is expressed in terms of the integration variables, namely: q° = (si + q 2 — S2)/2-Jsi. The functions 
Wi(si, S2,q 2 ) have the form (for more details see Appendix) 



q 2 4 

q 2 wi + —w 2 = -C(ml,ml 7 q 2 ), w x = 4(m x + m\ - q 2 ) - 16/3, 



w 3 



si(l - ai - a 2 ), 



(59) 
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with /3, a±, a-i given by eqs (36-38) in (TT) and C(mf, m 2 ,, q 2 ) defined in (|l7|). 

The electron spectrum dT/dE is obtained by integrating (|5^) over q 2 . Fig. || plots the results of the calculations. 
Here the effects of the subleading 1 /tuq orders are more pronounced but nevertheless lead only to a moderate change 
of the quark model box-diagram result. We also compare the quark model prediction with the electron spectrum in 
the free-quark decay process. 
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FIG. 6. Distribution dT/dE in B — » X c lt>t decays vs the lepton energy E: (a) ratio of the bound-to- free quark decay 
R(E) = (dT /dE)/(dTo/dE), lines same as in Fig. |, (b) Our QM result (solid) vs free-quark result (dot-dashed). 



VI. CONCLUSION 



We have proposed a new approach to the description of the quark-binding effects in the inclusive decays of heavy 
mesons. Our approach is based on the dispersion formulation of the relativistic quark model and allows one to express 
kinematical distributions in inclusive decays of heavy mesons in terms of the heavy meson soft wave function. This 
soft wave function describes long-distance effects both in exclusive and incluisve processes. 

Our main results are as follows: 

1. we have analysed the hadronic tensor in the quark model. We argue that the diagrammatic representation of 
the hadronic tensor in the quark model yields an expansion in inverse powers of the heavy quark mass as well 
as the standard OPE. However in distinction to the standard OPE which is based on the free quark decay as 
the LO process, the LO process in the quark model is described by the box diagram with a free Qq in the final 
state. This yields some specific features of the hadronic tensor calculated within the QM, both negative and 
positive. On the one hand, a consideration based on the box diagram alone reproduces the correct LO result 
but contains also 1/tuq corrections. These 1/tuq terms are cancelled against contributions of other graphs thus 
leading to the agreement with the standard OPE result. Hence the effects of the subleading diagrams should 
be taken into account for obtaining a consistent approach. On the other hand, the hadronic tensor calculated 
from the box diagram is a regular function in the kinematically allowed region of all variables, as well as all 
the subleading order I/tuq terms. This feature makes the quark model calculation of the hadronic tensor very 
suitable for describing the differential distributions; 

2. we have constructed the double spectral representation of the hadronic tensor for the B — > X c £Pi decay in 
terms of the B meson soft wave function and the double spectral density of the box diagram, and analysed its 
1 / tuq expansion in the case of the heavy-to-heavy inclusive transition. The spectral representation is further 
modified in order to take into account essential effects of the other diagrams contributing in subleading orders. 
Namely, the subtraction term in this dispersion representation is determined such that the 1 /tuq correction in 
the integrated SL rate is absent in agreement with the OPE result. Furthermore, a phenomenological cut is 
introduced in the spectral representation to bring the size of the 1 /tuq terms in the differential g 2 -distribution 
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at q 2 = and in the integrated rate into full agreement with the OPE. Thus our representation of the hadronic 
tensor obeys the OPE predictions in the regions where the latter are expected to be valid; 

3. we have obtained numerical results on differential distributions in inclusive B — > X c £p£ decays using the B- meson 
wave function and other quark model parameters previously determined from the description of exclusive meson 
transition form factors within the dispersion approach. So, basically our predictions are parameter-free. We 
notice that modifications of the spectral representations which take into account the subleading 1/toq effects 
within the box-diagram representation, introduce some uncertainties in our results. However, they do not affect 
our predictions strongly, and the main features of the inclusive distributions are determined by the soft meson 
wave function. Moreover, the size of the subleading corrections is in perfect agreement with the OPE result 
for the integrated rate, and we expect to describe also these subleading effects in differential distributions in a 
proper quantitative way. 

The proposed approach can be applied to the inclusive B — > X u ^ s transitions. In particular it is especially suitable for 
the description of the photon line shape in B — > X s ~f decays. However, certain subtleties in heavy-to-light transitions 
compared with the heavy-to-heavy decays emerge. They are mostly connected with the fact that in hcavy-to-light 
transition the kinematically allowed (^-interval of the hadron SL decay is larger than that of the quark decay, while 
in case of the heavy-to-heavy transitions the situation is just opposite and the q 2 -region of the quark decay is larger. 
This feature requires a detailed analysis of the <7 2 -region near zero recoil in heavy-to-light inclusive decays. 

It is also worth noting that our approach takes into account only non-perturbative effects in inclusive decays of 
heavy mesons. Perturbative corrections have been ignored. So, for comparing our results with the experimental 
differential distributions perturbative corrections should be also included into consideration. 

We are going to address these issues in a separate work. 
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VII. APPENDIX: CALCULATION OF THE DOUBLE DIFFERENTIAL DISTRIBUTION 



Here some technical details of calculaitng the d 2 T/dq 2 dE are provided. We start from the expression (|2^) which 
gives the double discontinuity of the box diagram. The trace corresponding to the V — A current has the form 

Sp (75 (ma - h)l5 (m 2 + hhni 1 ~ 7b) (mi + £1)7,/ (1 - "/5){m 2 + k 2 )j = 2 (si - (m 2 - m 3 ) 2 ) • w^, (60) 

where w)^,, is the trace over the free-quark loop 

«V = S P (( TO 2 + feh^ 1 _ 75) ("*i + hhv(l - 75)) = 8[fci p fc 2 „ + fe 2M fci„ - g^ukxki + ienvapfaaQp]- (61) 
One finds the following useful relation 

- 3 (<7 2 5^ - = ^( fc i fc 2 9 2 + 2fcig k 2 q) = - [(ml - m 2 ) 2 + g 2 (m 2 + m 2 ) - 2q i ] = ^C(m\, ml, q 2 ). (62) 

The integration over dk\dk 2 . . . in eq (|2^) yields 

(a \ tt6I(. . .) ( Sl - (m 2 - m 3 ) 2 ) 

( A ^J = \l/2( m 2 m 2 „ 2 N «V(Pl,<?) (63) 

where w^ v is represented in terms of the 'dispersion momenta' p\ and q — p\ — p 2 (p 2 = s±, p\ — s 2 , g 2 = q 2 ) as 
follows 

Wfi V {j>i,q) = -g^uWi H w 2 + it llva p^=qpWz H — u; 4 +q fJ .q,yW 5 . (64) 

si y Si Y Si 

Notice that 

- 3 (<? 2 .9^ ^ ?/^) = I 2 ™! + \ t w 2, (65) 



where 



A 1 / 2 (si,s 2 ,g 2 ) o si + g 2 -s 2 
^ = ^~7F ' q = ' ( 66 ) 



Comparing (p2j) and (pq) one finds 



q 2 w 1 + lf W2 ^^C(mlm 2 ,q 2 ). (67) 



As the next step, performing the convolution with the leptonic tensor gives the double differential distribution 
d 2 T _ G 2 F 



\ V 2i\ 2 f, 2, ,s 1 -(m 2 -m 3 ) 2 1/2 2 2 s F ( s i) f 1 drj ( q 2 \ , max 2 .v 

12"^ I dsiv (si) — ^ — A (si,ffl2 ' m3) ^" li tT° > E+ w e ^ < 52 (9 » 

x (2g 2 u;i( Sl , a 2 , g 2 ) + [4£(g° - E) - q 2 }w 2 { Sl , s 2 , q 2 ) + 2q 2 (2E - q Q )w 3 ( Sl , s 2 ,q 2 )) . (68) 



In this formula s 2 is connected with 77 through eq (|32| ) , and qo is given in terms of the integration variables by eq (|6J^ 

The eq ( |68| ) is based on the box diagram and also includes modifications which tune the size of the 1 /tuq corrections 
as explained in the text. 



By virtue of the relations 



J dE8(q° > E + q 2 /4E) = \q\, (69) 
J EdE6(q° >E + q 2 /4E) = q°\q\/2, (70) 
J E 2 dE8(q° >E + q 2 /4E) = ±\q\(3q 2 + <f). (71) 
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one finds 



J dE{2q 2 Wl + [AE(q° -E)- q 2 }w 2 + 2q 2 (2E - q a )w 3 }9{q° > E + q 2 / ' AE) = 2|g](gV + ^fw 2 ). (72) 
Thus, integrating the double differential distribution (|68|) over E and using (|67|) gives dT/dq 2 in the form (E3) 



dT G 2 p\V 2 l\ 2 , 2 2 2n /" , 2/ n«1 " ( TO 2 ~ TO3) 2 ,1/2, 2 2n 

= q« 3 3 Cim^m^q ) / dsi^ (si) A (si,m 2 ,m 3 ) 

(m 2 +'m 3 ) 2 

xFW^I T -^A 1 / 2 (si, s 2 , g 2 )6» (s 2 < s 2 nax (q 2 )) . (73) 

It is clear that for the calculation of dT / dq 2 we do not need to know all Wi , but only their linear combination ( |67| ) . 
On the other hand, for calculating the electron spectrum dT/dE we need all the functions Wi. The latter read 

a 2 4 

2 1 i /-if 2 2 2\ 

9 Wl I W2 = 3 6 ( m i> m 2>9 )> 

Wl = 8fcifc 2 - 16/? = A(m\ + m| - q 2 ) - 16/3, 

W3 = &y/si(l - Oil - a 2 ), (74) 

where /3, eti, a 2 are given by eqs (36-38) of Jll|. 

It might be interesting to note that at q 2 ^ 0, and using the reference frame q + = 0, p±±_ = (q 2 = —q 2 ^) one 
obtains 

= -^kl-^^j, l-ai-a 2 = l-x 3 , (75) 

where 2:3 and k± are the (+) and (_L) components of the spectator quark momentum, respectively. In the heavy 
meson one finds 

/3~A 2 , xa-A/mg. (76) 

Although at q 2 > the interpretation of (3 and ot\ + a 2 in terms of k± and X3 is not straightforward, the estimates j76| ) 
remain valid also at q 2 > 0. This means that in B — > X c lv(. decays our Wj's differ only slightly from the corresponding 
free-quark expressions. 
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